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Abstract. We have developed a new method for determining microscopically the five- 
dimensional quadrupole collective Hamiltonian, on the basis of the adiabatic self-consistent 
collective coordinate method. This method consists of the constrained Hartree-Fock-Bogoliubov 
(HFB) equation and the local QRPA (LQRPA) equations, which are an extension of the usual 
QRPA (quasiparticle random phase approximation) to non-HFB-equilibrium points, on top of 
the CHFB states. One of the advantages of our method is that the inertial functions calculated 
with this method contain the contributions of the time-odd components of the mean field, which 
are ignored in the widely-used cranking formula. We illustrate usefulness of our method by 
applying to oblate-prolate shape coexistence in ^^Kr and shape phase transition in neutron-rich 
Cr isotopes around N — 40. 



1. Introduction 

The five-dimensional (5D) quadrupole collective Hamiltonian is a powerful tool to describe 
large-amplitude deformation dynamics including triaxial shape degree of freedom as well as 
axially symmetric deformation. The 5D quadrupole collective Hamiltonian is characterized by 
seven functions: the collective potential, three vibrational and three rotational inertial functions 
(collective masses). To calculate the inertial functions, the Inglis-Belyaev (IB) cranking formula 
[1, 2] has been widely used. However, it is well known that IB cranking masses do not contain the 
contributions from the time-odd components of the mean field and underestimate the collective 
masses. 

In this presentation, we introduce a new method of determining microscopically the collective 
potential and the inertial functions which can overcome the shortcoming of the IB cranking 
formula. This method is derived from the adiabatic self-consistent collective coordinate 
(ASCC) method [3, 4], which is based on the adiabatic time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory and, with which one can extract the collective subspace (collective manifold) 
from the large-dimensional TDHFB configuration space. The method we employ in this study is 
the version of the ASCC method simplified by assuming the one-to-one correspondence between 
the collective manifold and the quadrupole deformation parameter space (/3,7). The main 
concept of this method is the local normal modes on top of the constrained Hartree-Fock- 
Bogoliubov (CHFB) states at every point of the (/3, 7) plane. In this method, we first solve the 
CHFB equations imposing the constraints on the proton and neutron numbers and deformation 



parameters (/3,7). Then, we solve local QRPA (quasiparticle random-phase approximation) 
equations, which are an extension of the usual QRPA to the non-HFB-equilibrium points. 
Hereafter, we cah the method "the constrained HFB plus local QRPA (CHFB+LQRPA) 
method." 

We shall show some results of the applications of the CHFB+LQRPA method to the oblate- 
prolatc shape coexistence in ^^Kr and to the shape phase transition in neutron-rich Cr isotopes 
around = 40. 



2. Theoretical framework 

Here, we briefly explain the theoretical framework of the CHFB+LQRPA method (see Ref. [5] 
for details). The 5D quadrupole collective Hamiltonian is written in terms of the magnitude 
(3, the degree of triaxiality 7 of quadrupole deformation, three Euler angles, and their time 
derivatives as 

^coU =7;ib + rrot + F(/3,7), (1) 

Tvib = ^L»^^(/3,7)/3' + i^/37(/5, 7)^7 + ^^77(/5, 7)7', (2) 
1 ^ 

T,ot = -Y,JkW,l)^l (3) 
k=l 

where Tyi^, T^ot and V represent the vibrational, rotational and collective potential energies, 
respectively. We determine the three vibrational inertial masses, three rotational moments of 
inertia, and collective potential in the collective Hamiltonian, by solving the CHFB+LQRPA 
equations. 

In the CHFB+LQRPA method, we first solve the CHFB equation to determine the collective 
potential V{P,j): 

<5(</.(/3,7)|i?CHFB(/3,7)l0(/5,7)) =0, (4) 

where i^cHFB = H -Y.r A(^^iV(^) - Ai^'"^^2m- Then, we solve the LQRPA equations on top 
of the CHFB states obtained above, 

5 (</.(/?, 7)1 [^chfb(/3,7),Q'(^,7)] - m,l)) = 0, (5) 

5 [i?CHFB(/3,7), \P0,1)] - Ci{P,j)Q\p,j) |</.(/3,7)) =0, (i = 1,2). (6) 

The vibrational inertial functions are calculated by transforming two LQRPA modes to the 7) 

degrees of freedom. We also solve the LQRPA equations for rotation to determine the moments 
of inertia. In this work, we adopt a version of the pairing-plus-quadrupole (P+Q) interaction 
including the quadrupole pairing interaction as well as the monopole pairing interaction and 
take two major harmonic oscillator shells as a model space. 

Finally, we solve the collective Schrodinger equation for the 5D quadrupole collective 
Hamiltonian quantized according to Pauli's prescription to obtain the excitation energies and 
collective wave functions, with which the electric quadrupole transitions and moments are 
calculated. 



3. Application to oblate-prolate shape coexistence in Kr 

In this section, we briefly show the numerical results of the application of the CHFB+LQRPA 
method to oblate-prolate shape coexistence in ^^Kr. The effective charges are adjusted to 
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Figure 1. Collective potential energy surfaces F(/3,7) in units of MeV, LQRPA vibrational 
inertial mass D^^(/3,7) and LQRPA rotational moment of inertia J'i{P,j) in unit of MeV"-*^ for 



(4ff '4ff ) = (0.658,1.658) such that the calculated result of B{E2;2f 0+) reproduces the 
experimental value for ^^Kr [7]. The detailed results of this calculation are shown in Ref. [6] 
together with the results for ^'''^^Kr. 

Figure 1 shows the collective potential y(/3, 7), the vibrational inertial function Z)^^(/3,7), 
and the rotational moment of inertia J7i(^, 7) for ^^Kr. The collective potential has two local 
minima: the oblate minimum is lower than the prolate one. The spherical shape is a local 
maximum. One can see that Z)/3^ indicates strong /? — 7 dependence and that the moment of 
inertia J7i deviates strongly from the irrotational moment of inertia. In Ref. [6] , we have seen the 
collective wave function is localized on the oblate (prolate) side in the 7) plane in the ground 
(excited) band and that the localization of the wave function develops with increasing angular 
momcntTim. The development of the localization is strongly related to the /3 — 'j dependence of 
the moments of inertia seen in Fig. 1(c). 

In Fig. 2, the excitation energies and B{E2) values for ^^Kr are shown together with the 
experimental data. We show the result obtained using the IB cranking masses for comparison. 
The excitation energies obtained with the LQRPA masses are lower than those obtained with 
the IB cranking masses and agree better with the experimental data except for the 2^ state. In 
particular, the observed O2" excitation energy which is close to the 2f excitation energy is well 
reproduced with the LQRPA masses. One can sec from the E2 transition strengths in Fig. 2 
that the shape-coexistence-like character becomes stronger with increasing angular momentum: 
the interband transitions between the initial and final states having equal angular momentum 
become weaker and weaker, which reflects the development of the localization of the vibrational 
wave function mentioned above. 

4. Application to shape transition in neutron-rich Cr isotopes 

In this section, we shall show the results of the application to the CHFB+LQRPA method to 
shape phase transition in neutron-rich Cr isotopes around A'^ = 40, where the development of 
deformation has been suggested by recent experimental data. More details will be discussed in 
our forthcoming paper [11]. 

In this calculation, we have determined the interaction parameters as follows: the parameters 
for ^^Cr are determined such that the calculated monopole pairing gaps and deformation at 
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Figure 2. Excitation spectra and B{E2) values calculated for '^Kr by means of the 
CHFB+LQRPA method (denoted LQRPA) and experimental data [8, 9, 10]. For comparison, 

results calculated using the IB cranking masses (denoted Inglis-Belyaev) are also shown. Only 
B[E2ys larger than 1 Weisskopf unit are shown in units of e^fm^. 



the HFB equilibrium point as well as the pairing gaps at the spherical shape reproduce those 
obtained with the Skyrme-HFB calculation with the SkM* functional using the HFBTHO code 
[12]. For the other nuclei 58,60,64,66qj.^ assume the simple mass number dependence according 
to Baranger and Kumar [13]. Wc take two major harmonic oscillator shells with iV = 3,4 and 
= 2, 3 for neutrons and protons, respectively. The single particle energies are determined 
from those obtained with the constrained Skyrme-HFB calculation at spherical shape. We scale 
them according to the effective mass of the SkM* functional m* /m = 0.79. For the calculation 
of the E2 transitions and moments, we have used a standard value of the effective charges 

(4ff^4S) = (0-5,l-5). 

We show the collective potentials for ^ Cr in Fig. 3. The location of the absolute minimum 

is indicated by the triangle. In ^^Cr, the absolute minimum is located at a nearly spherical 
shape. Although the minimum shifts to larger /3 in ^°Cr, the potential is extremely soft in the 
P direction. In ^^Cr, a more definite local minimum appears and the minimum gets still deeper 
in ^^Cr. The potential energy surfaces seem to indicate a shape transition from a spherical to a 
prolately deformed shape along the isotopic chain toward N = 40. 

We show in Fig. 4 the calculated excitation energies of the 2f and 4^ states, their ratios 
i?4/2 = E{4f)/E(2'^), B{E2;2'^ — )■ 0^*"), and the spectroscopic quadrupole moments of the 
2f states Q{2'^) calculated for 58-66(^j. comparison with the available experimental data. 
Both of the experimental excitation energies and R4/2 are reproduced well. The decrease in 
the excitation energies and the increase in R4/2 reflect the development of deformation with 
increasing the neutron number. The E2 transition and moments also indicate the enhancement 
of collectivity: the magnitude of B{E2) and Q(2^) increase as the neutron number increases and 
they take a maximum at A = 40. (Note that (5(2]'")'s are negative indicating prolate shapes.) 
To sum up, our result suggests the development of prolate deformation from AT = 34 to AT = 40 
and the largest collectivity at A" = 40. 





Figure 3. Collective potential energy surfaces 7) for ^^Cr in units of MeV. The triangle 
denotes the position of the absolute minimum. 



5. Concluding Remarks 

We have proposed the CHFB+LQRPA method for determining the inertial functions in the 5D 
quadrupole collective Hamiltonian, with which one can take into account the contributions from 
the time-odd components of the mean field to the inertial functions. We applied this method 
to the oblate-prolatc shape coexistence in the low-lying states of '^^Kr and shape transition in 
neutron-rich Cr isotopes around N = 40. The calculated results are in good agreement with the 
available experimental data. 

The CHFB+LQRPA method is based on the ASCC method and it can be used in conjunction 
with any interaction in principle. Nevertheless, in this study, we have employed a rather simple 
interaction, the P-l-Q model including the quadrupole pairing interaction, for simplicity. The 
implementation with modern energy density functionals, such as Skyrme energy functionals and 
density-dependent pairing interaction is an issue for future under progress. 
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Figure 4. Excitation energies of tlie 2'^ and states, tlieir ratios B{E2;2f — t- Oj^), 

and the spectroscopic quadrupole moments of the 2f states in comparison with the available 
experimental data [14, 15, 16, 17]. The B{E2) values are shown in Weisskopf units. 
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